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Abstract. We investigated the importance of spectral range 
and angular resolution for aerosol retrieval from multiangle 
photopolarimetric measurements over land. For this purpose, 
we use an extensive set of simulated measurements for differ- 
ent spectral ranges and angular resolutions and subsets of real 
measurements of the airborne Research Scanning Polarime- 
ter (RSP) carried out during the PODEX and SEAC 4 RS 
campaigns over the continental USA. Aerosol retrievals per- 
formed from RSP measurements show good agreement with 
ground-based AERONET measurements for aerosol optical 
depth (AOD), single scattering albedo (SSA) and refractive 
index. Furthermore, we found that inclusion of shortwave in- 
frared bands (1590 and/or 2250nm) significantly improves 
the retrieval of AOD, SSA and coarse mode microphysical 
properties. However, accuracies of the retrieved aerosol prop- 
erties do not improve significantly when more than five view- 
ing angles are used in the retrieval. 


1 Introduction 

The radiative impact caused by aerosols is widely consid- 

ered to be one of the largest uncertainties in the radiative 

forcing of the Earth climate (IPCC, 2013). Aerosols influ- 
ence the radiative balance of the Earth by scattering and 
absorbing solar radiation. A major issue for the evaluation 


of aerosols in climate models is the limited accuracy and 
detail of global measurements of the aerosol properties re- 
quired for climate modeling, such as particle size distribu- 
tion, shape, complex refractive index and columnar concen- 
tration. Orbital remote sensing instruments that perform mul- 
tiangle, multispectral photopolarimetric measurements have 
the capability to provide these aerosol properties with the re- 
quired accuracy (Mishchenko and Travis, 1997; Hasekamp 
and Landgraf, 2007). The advantages of polarimetric mea- 
surements over intensity-only measurements come from the 
sensitivity of polarization of light and its spectral and angular 
dependence to particle size, shape and refractive index. The 
Polarization and Directionality of the Earth’s Reflectance 
(POLDER), the most recent in-orbit polarization-measuring 
satellite instrument, has demonstrated the importance of po- 
larization measurements for the retrieval of aerosol proper- 
ties (Hasekamp et al., 2011; Dubovik et ah, 2011). Limita- 
tions of the POLDER instrument are its inherent limited po- 
larimetric accuracy and its relatively small spectral range. To 
overcome these limitations, the Aerosol Polarimetry Sensor 
(APS) was scheduled for launch on the Glory mission. This 
instrument was designed to view a scene at ~ 250 viewing 
angles in nine spectral bands ranging from 410 to 2250 nm. 
Due to the tragedy of the Glory launch failure and the de- 
commissioning of POLDER, space-borne multiangle photo 
polarimetric measurements are currently unavailable. How- 
ever, the Research Scanning Polarimeter (RSP) that has sim- 
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ilar characteristics as APS continues to perform airborne 
measurements in many field campaigns. The collected data 
will facilitate algorithm development and allow us to derive 
geophysically relevant aerosol and cloud properties (Waquet 
et al., 2009a). 

RSP provides measurements of a ground scene at ~ 152 
different viewing angles (±60° from nadir) in nine spec- 
tral bands from visible/near-infrared bands to shortwave in- 
frared bands (410, 470, 550, 670, 865, 960, 1590, 1880 and 
2250 nm). The measurements provide the first three Stokes 
parameters, /, Q and U , in the nine spectral channels with 
a radiometric uncertainty of ~ 2.0% and a polarimetric un- 
certainty of ~0.5% (Cairns et al., 1999). Given its large 
spectral range and number of viewing angles, RSP measure- 
ments are well suited to study the importance of spectral 
range and angular resolution for the retrieval of aerosol prop- 
erties. These are two important design parameters for future 
polarimetric instrumentation. 

In this paper we apply an extended version of the SRON 
aerosol retrieval algorithm (Hasekamp et al., 2011) to sev- 
eral measurement subsets of RSP measurements over land 
in order to investigate the importance of spectral range and 
angular resolution of multiangle photopolarimetric instru- 
mentation for aerosol retrieval. In addition, we also study 
the above-mentioned aspects using an extensive set of syn- 
thetic measurements. Although studies using synthetic multi- 
angle photopolarimetric measurements have been published 
before (Hasekamp and Landgraf, 2007; Hasekamp, 2010; 
Knobelspiesse et al., 2012; Ottaviani et al., 2013), they were 
all based on linear error propagation and a limited number 
of aerosol scenarios. Instead, here we perform a full iterative 
retrieval on the synthetic measurements and consider a large 
number of aerosol scenarios. Therefore, the synthetic study 
included in this paper is an extension of the synthetic studies 
mentioned above. 

In order to exploit the full spectral range of RSP, we extend 
the SRON aerosol retrieval algorithm to be capable to cope 
with spectrally dependent refractive indices. The paper is or- 
ganized as follows: Sect. 2 describes the data used in the pa- 
per, Sect. 3 describes the retrieval method. Sect. 4 shows the 
results based on synthetic measurements and Sect. 5 shows 
the results based on RSP measurements. Finally, Sect. 6 con- 
cludes the paper. 


2 Data 

The RSP data used in the paper are obtained during the 

PODEX (Polarimeter Definition Experiment) campaign and 
the SEAC 4 RS (Studies of Emissions and Atmospheric Com- 
position, Clouds and Climate Coupling by Regional Sur- 
veys) campaign in 2013 (Table 1). While the objectives for 
both campaigns include studies for future instrument design, 
SEAC 4 RS also seeks to improve our understanding of feed- 
backs of aerosols on meteorology and climate (Toon, 2013). 
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Figure 1. Measurements of individual ground pixel and effective 
5 km averaged pixel at 470 nm. Oscillations in angular measure- 
ments are due to imperfect co-registration of different viewing an- 
gles. 

The PODEX and SEAC 2 * 4 RS measurements were acquired 
over the continental USA with RSP instrument mounted on 
an ER-2 high-altitude aircraft which flies at an altitude of 
approximately 20km with a speed of around 200ms _1 . Dur- 
ing the SEAC 4 RS campaign, a number of AERONET sites 
co-located with the ER-2 flight track are available; thus we 
can compare aerosol properties obtained from RSP measure- 
ments with corresponding AERONET aerosol products. 

For both campaigns, we focus on measurements of the 
RSP bands at 410, 470, 550, 670, 865, 1590 and 2250nm. 
The 960 and 1880nm bands are excluded, because their main 
purposes are water vapor estimation and cirrus cloud screen- 
ing. For each spectral band, we only use viewing angles in 
the range [—60, 40°] because some of the viewing angles are 
blocked by the ER-2 aircraft. We use the data that have al- 
ready been aggregated to ground level. Due to imperfect co- 
registration of different viewing angles, the angular measure- 
ments show oscillations larger than the measurement error 
for inhomogeneous scenes. In order to reduce this effect we 
average RSP measurements over a distance of 5 km so that 
misregistrations between viewing angles become small com- 
pared to the effective pixel size. Figure 1 shows an example 
of a RSP measurement for an individual ground pixel and for 
an effective 5 km averaged pixel. 

Here, information on temperature, pressure and relative 
humidity profiles in atmosphere is interpolated to a partic- 
ular date, time and location of the RSP measurements using 
NCEP reanalysis data (Kistler et al., 2001). 
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Table 1. Information of used RSP measurements 


Case 

Campaign 

Co-located AERONET site 

Longitude and latitude 

Date 

1 

PODEX 

Fresno_2 

119.77°W, 36.78°N 

20130116 

2 

PODEX 

Fresno_2 

119.77°W, 36.78°N 

20130120 

3 

PODEX 

DRAGON_Porterville 

119.05°W, 36.03°N 

20130120 

4 

PODEX 

DRAGON_Drummond 

119.74°W, 36.70°N 

20130120 

5 

PODEX 

DRAGON_ Visalia 

119.39°W, 36.31°N 

20130122 

6 

SEAC 3 4 RS 

Calipso_Carthage 

94.06°W, 32.06°N 

20130909 

7 

seac 4 rs 

Upper_Buffalo 

93.20°W, 35.82°N 

20130906 

8 

seac 4 rs 

Upper_Buffalo 

93.20°W, 35.82°N 

20130827 

9 

seac 4 rs 

Leland_HS 

90.89°W, 33.40°N 

20130909 

10 

seac 4 rs 

Leland_HS 

90.89°W, 33.40°N 

20130830 

11-14 

seac 4 rs 

Baskin 

91.73°W, 32.28°N 

20130909 



1000 1500 2000 2500 


Wavelength (nm) Wavelength ( nm ) 

Figure 2. Imaginary (left panel) and real parts (right panel) of the refractive index as a function of wavelength. Data taken from D’ Almeida 
et al. (1991). 



3 Methods 

3.1 Forward model 

In the PODEX and SEAC 4 RS campaigns, the airborne RSP 
measurements are performed at an altitude of about 20 km. 
To model the radiance and state of polarization of light of 
a given wavelength measured at the RSP flight altitude, we 
consider a macroscopically isotropic, plane -parallel multi- 
layered atmosphere bounded below by a reflecting surface. 
Each atmosphere layer is characterized by scattering and ab- 
sorption optical thickness and the scattering phase matrix P. 

For RSP shortwave infrared bands, gaseous absorption 
need to be taken into account. The absorption in RSP 1590 
and 2250 nm bands can be represented by a model that is 
derived using k distributions for each 1 nm that are then inte- 
grated with a solar weighting over the RSP bandpasses. Here 
the total absorption optical depth depends on the total column 
optical depth of well-mixed gases and water vapor. It should 


be noted that water vapor absorption is relatively small when 
compared with that of the well-mixed gases CO 2 at 1590 nm 
and methane at 2250 nm. In terms of profile distributions, the 
integrated optical depths from the top of the atmosphere to 
a pressure level p for the well-mixed gases and water vapor 
are proportional to (p/ po) 2 and ( p/ po) 4 , respectively, where 
po is the surface pressure. 

The optical properties of aerosols depend on the size, 
shape and type of aerosols. There is sufficient evidence that 
mixtures of spheroids allow rather accurate fitting of mea- 
sured spectral and angular dependencies of observed inten- 
sity and polarization (Dubovik et al., 2006). Here we con- 
sider scattering properties of both spherical and spheroidal 
aerosol particles. The optical properties of spherical and 
spheroidal aerosol particles are calculated using the tabulated 
kernels of Dubovik et al. (2006). Here we describe particles 
as a mixture of spheres and spheroids with the aspect ratio 
distribution proposed by Dubovik et al. (2006). 
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Figure 3. Comparison between accounting for (top panel, for / 2 < 2.0, 35.3% cases converged) and not accounting for (bottom panel, 
75.1 % converged) spectral variation of refractive index in synthetic retrievals. The correlation coefficients (cor), mean error (Am) and mean 
absolute error (< Am >) between retrieved and true values of refractive index are included. 


For the characterization of surface reflection, the Rahman- 
Pinty-Verstraete model (RPV) (Rahman et al., 1993) is used 
with three independent parameters: po- g and k (see Ap- 
pendix). The model successfully accounts for the observed 
variability of reflectance measurements in laboratory and 
field conditions, ranging from bare soil to full canopy cover, 
in both the visible and the near-infrared bands (Rahman 
et al., 1993). To account for surface polarized reflectance, 
we use a modified Fresnel model introduced by Litvinov 
et al. (2011) (see Appendix). It has been shown that the Bidi- 
rectional Reflectance and Polarization Distribution Function 
(BRDF and BPDF) used here can fit the RSP measurements 
with sufficient accuracy (Litvinov et al., 2011). 

We solve the vector radiative transfer equation in the Earth 
atmosphere using the method as described by Hasekamp and 
Landgraf (2002, 2005a). To solve the integration over az- 
imuth angle, the radiative transfer equation is decomposed 
into a series of Fourier components. To reduce the number of 
terms needed in the Fourier series, the intensity vector orig- 
inating from single scattering is calculated analytically. To 
efficiently treat highly peaked phase function we apply the 
“MS method” by Nakajima and Tanaka (1988). Then, an it- 
erative Gauss-Seidel method is applied to solve the radiative 
transfer equations (Herman and Browning, 1965; Landgraf 
et al., 2001; Schepers et al., 2014). 


3.2 Retrieval method 

The aim of aerosol retrieval is to convert observed measure- 
ments into information about physical properties of aerosols, 
i.e., 

y=F(x) + e y . (1) 

Here y is the measurement vector that contains the multiple- 
wavelength multiple -viewing angle measurements of radi- 
ance and degree of linear polarization (DOLP), F is the for- 
ward radiative transfer model and e v is an error term. To 
define the state vector x, we employ the bimodal lognor- 
mal function to describe the atmospheric aerosol size dis- 
tributions including a fine and coarse mode indicated by su- 
perscripts ,v and c. Here, each mode is characterized by the 
effective radius r e rr, the effective variance i> e ff and aerosol- 
loading N that is the number of aerosol particles per square 
meter. Spherical and spheroidal aerosol particles are distin- 
guished by fraction ratio p, which indicates the percentage 
of spheroid particles. We consider a Gaussian height distri- 
bution h(z) — Aexp(— 41n2(z — z c ) 2 /z w ) with center height 
Z c and width z,„ for layer at height z and include z c as one 
additional unknown parameter in our retrieval. 

Additionally, the complex refractive index m — m r + im, , 
is needed to characterize aerosols of each mode. Figure 2 
shows the spectral dependence of the most important aerosol 
types (D’Almeida et al., 1991). It can be seen in Fig. 2 that 
for the aerosol types that predominantly occur in the fine 
mode (biomass burning, industrial, soot, sulfate and water 
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number of viewing angles 


number of viewing angles 


Figure 4. The results of mean absolute errors of aerosol optical depth ( 1550 ), single scattering albedo (&> 55 o), aerosol size parameters (V^g, 
u eff’ r eff alM ^ ^eff) an< ^ re f ract i ve index (m®, m®, m £ and m?) as a function of the number of viewing angles between retrieved and true 
values for four different RSP measurement configurations. Soil and vegetation surfaces are used in simulation. Aerosol types in Fig. 2 are 
all included in simulation and spectral variation of refractive index are included in both simulation and retrieval. The shaded areas indicate 
the target accuracy on aerosol parameters formulated by Mishchenko et al. (2004). For the target accuracy of AOD, we use 10 % of the mean 
AOD, which is 0.68. 
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soluble) the spectral dependence of the real part of the re- 
fractive index is relatively flat in the spectral range where 
they are most optically active (< 1200nm). Although the 
imaginary part of the refractive index for the fine mode 
types shows a spectral dependence for some types, the val- 
ues for these types are generally small, so that such wave- 
length dependences will not significantly affect optical mea- 
surements in the visible spectral range. This may be differ- 
ent when “brown carbon”-type aerosols are included in the 
analysis (Feng et al., 2013). For the coarse mode the situa- 
tion is different, since for this mode the real part of the re- 
fractive index generally shows a spectral dependence in the 
range >~ 1200nm where coarse mode aerosols are optically 
active. Furthermore, for dust the imaginary part of the re- 
fractive index increases strongly in the blue/UV part of the 
spectrum. Here, the spectral refractive index of coarse mode 
aerosols are treated as a mixing of three classes refractive in- 
dex spanning the wavelength range in Fig. 2 and thus can be 
described by a linear combination of them. If so, we have six 
parameters to be fit for both the real and imaginary part of 
refractive index of coarse mode aerosols. However, we find 
that these six parameters are strongly intercorrelated, so in 
fact less than six parameters are needed. The well-known 
principal component analysis method can be used to deter- 
mine the most important spectral signatures in the refractive 
index (Smith, 2002). For the real part of complex refractive 
index, at xn matrix M with n > t is constructed as 

M = [mi, m 2 , ■ • •, (2) 

where T indicates the transposed vector and row vectors are 
the mean-subtracted real part of the refractive indices at n 
different wavelengths for t types of representative coarse 
mode aerosols. From the covariance matrix of M, eigenvec- 
tors m A can be derived in order of significance indicated 
by eigenvalues. We order the vector m c j according to their 
eigenvalues and map them into the original data space by 
multiplying with M. Consequently, the original row vec- 
tors in matrix M or their linear combinations can be ex- 
actly reconstructed by wavelength-dependent vector m‘{ and 
wavelength-independent value /h 111(1 through 

t' 

m — m md + «,■ • m c j , (3) 

i = 1 

where t' = t and a; is the weight of vector m c ] . For coarse 
mode aerosols in Fig. 2, we found that for t' = 1 we can suf- 
ficiently describe the spectral features in the refractive index 
both for the real and the imaginary part. 

In total, for coarse mode aerosols there are eight state vec- 
tor elements of which two correspond to size (r e ff and n e ff), 
four correspond to refractive index (mj. , a r , m™ d and a,) 
and two others correspond to aerosol-loading N and spheroid 
fraction p, respectively. The fine mode aerosols are described 
in the same way with five state vector elements, i.e., without 
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the parameters describing spectral dependence of the refrac- 
tive index and the spheroid fraction. In the retrieval, we also 
retrieve the central height z c of the aerosol Gaussian height 
distribution. As a consequence, the state vector x contains 14 
aerosol-related elements in total. In addition to these aerosol 
parameters we also include the surface reflection parameters 
Pq, g and k of RPV model and surface polarized reflectance 
parameters a and a of the modified Fresnel model in the state 
vector. 

To retrieve the state vector x from RSP measurements, 
an inversion algorithm based on the Phillips-Tikhonov reg- 
ularization is employed (Hasekamp and Landgraf, 2005b) 
which has been successfully applied to aerosol retrieval from 
PARASOL measurements (Hasekamp et al., 2011). The in- 
version algorithm chooses the solution x that solves the min- 
imization problem 

x = minflISy 2 (F(x) - y) || 2 + y II W(x - x a ) || 2 ), (4) 

X 

where jc a is a priori state vector, W is a weighting matrix 
and S v is the measurement error covariance matrix that con- 
tains the measurement error estimate. The weighting matrix 
W is chosen diagonal such that W,j — \/x a j, making the 
side constraint dimensionless. The regularization parameter 
y is found from a method similar to the L curve (Hansen 
and O’Leary, 1993), in which we evaluate the least squares 
norm for different trial values of y using an approximate for- 
ward model and choose y for which the least squares norm 
is smallest. 

To obtain an appropriate first guess state vector for the 
iterative retrieval procedure, we perform a retrieval using a 
lookup table based on 48 aerosol modes for the fine mode 
and 24 for the coarse mode. Table 2 shows the parameters of 
each mode. 

For surface reflectance, the parameters of the RPV model 
used for generating the lookup table at each wavelength are 
shown in Table 3. For surface polarized reflectance, we use 
one mode of the modified Fresnel model in Table 3. The com- 
bination of fine and coarse mode and surface parameters that 
yields the best fit to the measurement is used as first guess 
and a priori in the full retrieval procedure. 

4 Aerosol retrieval 
4.1 Synthetic retrieval 

First, we perform our investigations on a set of 1000 syn- 
thetic measurements in the principal plane at solar zenith an- 
gle 60.0°. Each synthetic measurement includes 49 equally 
sampled viewing angles (±60.0° from nadir) and seven of 
the nine RSP wavelength bands as described in Sect. 2. The 
relative azimuthal angles are 0.0 and 180.0° for positive and 
negative viewing zenith angles, respectively. Here we use 
the convention that an azimuth angle of 180.0° corresponds 
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Figure 5. Comparison between retrieved and true values of aerosol optical depth ( 7550 ), single scattering albedo ( 01550 ) and aerosol pa- 
rameters (r®g-, m f’ m i> r eff’ w eff’ nif and mf) when the viewing angle number N = 5 using the full spectral range 410-2250nm. The 
correlation coefficients (cor), mean error (A) and mean absolute error (< A >) are included. 


Table 2. Parameters of fine and coarse modes used in lookup table retrieval. 


Modes 

refflpm] 

^eff 

m Y 


P T 

Fine mode 

0.05,0.10, 0.20, 0.30 

0.19 

1.37. 1.45, 1.52, 1.60 

0.001,0.020,0.100 

Q j Q 0.01, 0.15, 0.25, 0.5, 0.8, 1.0, 1.5, 3.0, 5.0 

Coarse mode 

1.00,2.50 

0.43 

1.3S, 1.44, 1.50 

0.001, 0.050 



80 100 120 140 160 180 

Scattering angle [°] 


Figure 6. One example of RSP measurements and the correspond- 
ing fit results for DOLP in retrieval. Triangle symbols denote the 
RSP measurements and solid lines with plus are fits using the for- 
ward model in Sect. 3. The date, AERONET site (longitude, lat- 
itude), number of viewing angles, AOD, mean angular deviation 
from the principal plane and solar zenith angle are 9 September 
2013, Baskin (91.73°W, 32.28°N), 60, 0.27 (440 nm), 1.99° and 
66.5°, respectively. 


to backward scattering. To get measurement configurations 
with different angular resolution, viewing angles are evenly 
sampled from —60.0 to 60.0° according to the resolution ex- 
cept in the measurement configurations with two viewing an- 
gles for which we use nadir and —60.0°. Aerosol size distri- 
bution parameters used here are randomly distributed in their 
acceptable range for both fine mode (0.1 < < 0.30 pm, 

0.1 < i>g ff < 0.30) and coarse mode (0.60 — r eff — 3 .50 pm, 
0.4 < iA |T < 0.6). To account for the spectral dependence, the 
refractive indices for each mode are linear combinations of 
the refractive index spectra as shown in Fig. 2. The AOD 
r for each mode varies randomly in the range 0.0 < 1550 < 
0.70. Aerosol layer center height are fixed to be 2000 m with 
a width of 2000 m. We consider soil and vegetation land sur- 
faces of RPV model and the modified Fresnel model with 
parameters presented in Table 4 (Litvinov et al., 2011). Pa- 
rameters po at 410, 470, 550 and 865 nm are obtained by lin- 
ear interpolation using soil and vegetation surface reflectance 
provided by ASTER spectral library (Baldridge et al., 2009). 

The forward model as described in Sect. 3 was used to 
simulate the intensity vector at the altitude of RSP (20 km) 
for all 1000 combinations of input parameters. In our simu- 
lations, we add an random noise of 2.0 % on intensity I and 
0.2 % on DOLP P to account for measurement uncertainties. 
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Table 3. Parameters of RPV model and the modified Fresnel model used in lookup table retrieval. 


A[nm] 

P0 

g and k 

a , a and k y 

410 

470 

Other bands 

0.02, 0.04, 0.06, 0.10, 0.15 
0.02, 0.05, 0.08,0.15, 0.25 
0.10, 0.15, 0.20, 0.30, 0.45 

g = —0.1, k = 0.2 or g = 0.5, k = 1.0 

a = 1.0, a — 2.0 and k y = 0.69 



Figure 7. Comparison between the retrieved AOD of RSP measurements and AERONET at 440, 675 and 870 nm. The correlation coefficients 
(cor) at each wavelength are included. 



AERONET AERONET 

Figure 8. Comparison between the retrieved fine and coarse 
modes effective radius (r£j. and of RSP measurements and 
AERONET. The correlation coefficients (cor) for each mode are in- 
cluded. Plus indicates cases where 7575 <0.1. 

4.1.1 Spectral dependent refractive index 

We first test the effect of spectral variation of the refractive 
index in the retrieval. A measurement configuration is used 
with a wavelength range of 41 0-2250 nm and 13 viewing an- 
gles. We consider retrievals to be successfully converged if 
the x 2 of the fit < 2.0 with x 2 defined as ^^7^_ 1 ( - v, '~ f ’' ) 2 , in 
which N is the number of measurements, y, is the measure- 
ment, Fj is the simulated result and n, is the error of the mea- 
surement. We found that this x 2 threshold yields a reasonable 
balance between accuracy and number of converged cases. 
Figure 3 compares retrieved refractive indices (at 670 nm) 
of fine and coarse modes with true values for retrievals us- 
ing either a constant or a wavelength-dependent refractive 
index. The general trend is that for fine mode aerosols, as 


expected, both approaches yield similar results; for coarse 
mode aerosols, considering the spectral variation improves 
the agreement between retrieved and true values for both the 
real and the imaginary part of the refractive index. Most im- 
portantly, the number of retrievals that converge with a suf- 
ficiently good fit to the measurement (x 2 < 2 . 0 ) more than 
doubles from 35.3 to 75.1 %. 

4.1.2 Comparison between measurement 
configurations 

To investigate the importance of spectral range and angu- 
lar resolution, we investigate retrieval errors on T 550 (AOD), 
CU 550 (single scattering albedo, SSA), r® ff , r£ ff , Ug ff , m x , 
m ?, m x and m'r using different measurement configurations. 
Here we focus on four different wavelength ranges, namely 
410-670, 410-865, 410-1590 and 410-2250nm. The num- 
ber of viewing angles changes between 2 and 49. Although 
the number of cases that converge to x 2 < 2.0 may change 
per measurement configuration, we slightly adjust the x 2 
threshold for each measurement configuration so that the 
convergence percentage is 75 % for all of them. 

Figure 4 shows the results of mean absolute errors of each 
parameter as a function of the number of viewing angles ( N ) 
for the four wavelength ranges. It is important to note that 
the results remain virtually unchanged when we use a dif- 
ferent realization of 1000 combinations of aerosol param- 
eters or when we increase the number of realizations. The 
shaded areas indicate the target accuracy on aerosol param- 
eters formulated by Mishchenko et al. (2004). First of all, 
for all four wavelength ranges the retrieval errors are signif- 
icantly decreased when the number of viewing angles N is 
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Table 4. Parameters of RPV model and the modified Fresnel model used in synthetic retrieval. 




RPV model 

Modified Fresnel model 

Surface types 

PO 

8 k a 

<7 

ky 


Soil 

0.001(410nm), 
0.014(550 nm), 
0. 106(865 nm), 
0.116(2250nm) 

0.007(470 nm), -0.097 0.746 

0.071 (670 nm), 

0. 159(1590 nm). 

4.260 

0.589 

0.700 

Vegetation 

0.013(410nm), 
0.024(550 nm), 
0. 120(865 nm), 
0.060(2250 nm) 
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Figure 9. Comparison between retrieved single scattering albedo (SSA) and composed refractive index of RSP measurements and those 
provided by AERONET at 675nm. (Assa), (A m t ) and (Am,-) are the mean absolute difference of SSA and real and imaginary part of the 
refractive index between RSP and AERONET. Plus indicates cases where T 675 <0.1. 



number of viewing angles 


Figure 10. Retrieved AOD mean absolute error dependence on 
number of viewing angles for different spectral ranges of RSP data. 


increased from 2 to 3. The decrease after N > 3 is moderate 
in compared with the change from 2 to 3. The mean errors do 
not improve significantly when the number N is greater than 
~ 5. Moreover, Fig. 4 also shows the advantage of the mea- 
surement configurations with spectral ranges 410-1590 and 
41 0-2250 nm, as they yield smaller retrieval errors than mea- 
surement configurations with a smaller spectral range. When 
compared with the 410-670 and 410-865 nm cases, the mean 


error on the AOD at 550 nm T550 is reduced from ~ 0.1 and 
~ 0.06 to ~ 0.04 as shown in the upper left figure. The mean 
error on the SSA at 550 nm W550 is less than 0.02, while for 
the 410-670nm case the errors are around 0.03. The mea- 
surements of 410-1590 and 41 0-2250 nm yield smaller er- 
rors especially for the coarse mode. For example, the mean 
error on the effective radius is reduced from ~ 0.3 to 
~ 0.2 pm and the mean error on the real and imaginary part 
of the refractive index is also reduced. This is attributable to 
the fact that shortwave infrared bands are sensitive to coarse 
mode aerosols. For fine mode aerosols, the measurements of 
410-865 nm perform almost equally well as retrievals includ- 
ing shortwave infrared bands. Finally, the 410-670nm case 
usually produces the largest mean error for all parameters. In 
contrast, the 410-2250nm case can meet the target accuracy 
requirements indicated by shaded area for nearly all parame- 
ters. When we compare retrievals in the range 410-1590nm 
with retrievals in the range 410-2250nm, we see that the ac- 
curacy of the aerosol optical depth retrievals is only mini- 
mally reduced with the latter setup, while for the other pa- 
rameters the performance is similar. 

Figure 5 shows scatterplots of retrieved vs. true values for 
T550, a>55o and parameters of both aerosol modes when the 
viewing angle number N = 5 using the full spectral range 
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41 0-2250 nm. The correlation coefficients between true and 
retrieved values for each parameter are 0.97 ( 7550 ), 0.98 
(tt> 55 o), 0.97 (fine mode effective radius r* ff ), 0.71 (fine mode 
effective variance u eff)> 0.85 (fine mode real part of the re- 
fractive index m*), 0.97 (fine mode imaginary part of the re- 
fractive index m?), 0.88 (coarse mode effective radius rS f ), 
0.21 (coarse mode effective variance 0.93 (coarse mode 
real part of the refractive index m c r ) and 0.38 (coarse mode 
imaginary part of the refractive index mf). The retrievals of 
the imaginary part of the coarse model aerosol can be poor 
when fine mode aerosols dominate. However, as shown in 
Fig. 5, the retrieved single scattering albedo values agree well 
with true values, which indicates good retrievals on total ab- 
sorption. 

4.2 Aerosol retrieval using the RSP data 

We applied our algorithm described above to RSP data over 
land close (distance < 20 km) to AERONET stations so that 
we can compare our results with independent measurements. 
The data were pre-processed as described in Sect. 2. Re- 
trievals are made only on cloud-free measurements because 
clouds are not included in the forward model. Time differ- 
ences between the RSP measurements and AERONET prod- 
ucts are restricted to less than 1 h. We use the level 2 di- 
rect sun products for AOD comparison and level 1.5 in- 
version products for refractive index and single scattering 
albedo comparison. In the AERONET products, aerosol op- 
tical depth values are provided at wavelengths from 340 to 
1640nm with an estimated uncertainty of 0.01-0.02, which 
are very well suited for validation of the AOD retrieved by 
RSP. 

For the aerosol properties from the AERONET inversion 
product it should be noted that the corresponding accuracies 
are less well quantified and expected to be similar or larger 
than those of the RSP retrievals. So, for SSA, refractive in- 
dex and effective radius, the results presented below should 
be considered as a comparison rather than a validation. Apart 
from this, AERONET retrieves one refractive index whereas 
we retrieve refractive indices for both the fine and coarse 
mode. To get a general impression on the agreement between 
AERONET and RSP retrieval, we compose a refractive index 
m comp that is a weighted average by AOD r of each mode: 

m s T s + m c T c 

^comp — : • (5) 

r s + r c 

The weighted refractive index m comp should reflect to some 
extent the retrieval capability on this quantity. 

In our RSP retrieval, we use five viewing angles which are 
evenly sampled within the available range of each measure- 
ment and consider cases that passed a goodness of fit crite- 
rion under X 2 < 10.0. The x 2 threshold here is larger than 
the one we used in synthetic retrieval, because some mea- 
surements are still hampered by angular oscillations caused 
by inhomogeneity of the underlying surface or atmosphere. 
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Besides this, we include only cases where the sampled scat- 
tering range includes the 85.0°-155.0° interval because re- 
trieval errors are considerably larger for smaller angular 
ranges (Hasekamp and Landgraf, 2007). Figure 6 shows one 
example of RSP measurements and the corresponding fit re- 
sults for DOLP. Generally, the fits are good for each spec- 
tral band. In some cases, the measurements can be found to 
fluctuate between some viewing angles, thus making the fits 
worse (see the fit for 670 nm in Fig. 6). 

Figure 7 shows a comparison between the retrieved 
AOD from RSP measurements and the AOD provided by 
AERONET at 440, 675 and 870 nm. Here we show re- 
trieval results using the wavelength range 410-1590nm. 
From Fig. 7 it follows that retrieved AODs at the three wave- 
lengths are in good agreement with the AERONET results. 
The correlation coefficients are 0.96 for 440 nm, 0.97 for 
675 nm and 0.92 for 870 nm. 

In Fig. 8 the retrieved effective radii of fine mode aerosol 
agree well with those of AERONET. The bias in coarse mode 
is mainly attributable to the fact that the retrieved AODs of 
the coarse mode are relatively small (mostly less than 0.1), 
which makes the retrieved values less reliable. 

Figure 9 shows the values of the SSA, refractive index 
m and the corresponding AERONET values at 675 nm. For 
the SSA, the retrieval results agree reasonably well with 
AERONET values. The mean absolute difference in SSA be- 
tween AERONET and RSP is 0.04. Under the assumption 
that the AERONET refractive index can be approximated by 
the weighted average of the two modes, the retrieved values 
are relatively close to those provided by AERONET. As ex- 
pected, the deviation between RSP and AERONET is larger 
for cases with small AOD (rg 75 < 0.1) as indicated by the 
color code in Figure 9. 

The retrieved center heights z c varied in the range of 0.97 
and 2.66km and for each case a similar height range is re- 
trieved by the Cloud Physics Lidar (CPL) (not shown). A 
more extensive comparison between retrieved heights from 
RSP and those from CPL will be subject to future investiga- 
tion. 

Figure 10 shows the difference in AOD between RSP and 
AERONET as a function of the number of viewing angles 
used in the retrieval. Due to the limited number of cases for 
which AERONET SSA and refractive index measurements 
are available, we restrict the analysis to AOD only. We should 
notice that with the increase of N the retrieval performance 
sometimes becomes even a bit worse, which is due to the fact 
that the measurements at some viewing directions can be dis- 
torted by the inhomogeneity of the land surface. In general, 
the results of Fig. 10 confirm our findings from the study 
with synthetic measurements: (1) the strongest improvement 
in accuracy is seen when N increases from 2 to 3; (2) the re- 
trieval performances are not significantly further improved if 
more than five viewing angles are used in the retrieval; (3) re- 
trievals perform significantly better when shortwave infrared 
bands are included. 
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5 Conclusions 

In this paper we performed aerosol retrievals over land from 
multiangle photopolarimetric measurements produced syn- 
thetically and obtained by the airborne RSP instrument. We 
adjusted the SRON aerosol retrieval algorithm to cope with 
spectrally dependent refractive indices so that it can be ap- 
plied to all RSP spectral channels. 

We investigated the capability of different retrieval setups 
with varying spectral ranges and angular resolutions by con- 
sidering 1000 synthetic measurements for different aerosol 
scenarios. We find that the errors on the retrieved parameters 
were significantly decreased when including multiple view- 
ing geometries. When more than five viewing angles are in- 
cluded in the retrieval, virtually no significant improvement 
was found anymore. It should be noted that for future instru- 
mentation a higher number of viewing angles is still desirable 
in order to separate scattering by aerosols and cloud droplets 
by means of the rainbow (Waquet et al., 2009b; Hasekamp, 
2010) and for the retrieval of cloud microphysical quanti- 
ties (Alexandrov et al., 2012). Moreover, we demonstrated 
the advantage of including 1590nm shortwave infrared bands 
in the retrieval over retrievals restricted to the visible spec- 
tral range. This improvement is in particular significant for 
coarse mode aerosol parameters such as effective radius and 
refractive index. Here, the spectral range 410-1590 and 410- 
2250 nm showed similar performance. 
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For retrievals from real RSP measurements, the retrieved 
AOD agrees well with AERONET products. The corre- 
lation coefficients are 0.96 (440nm), 0.97 (675 nm) and 
0.92 (870 nm). The retrieved SSA also agrees well with 
AERONET although it should be noted that the number of 
comparison points is limited. The retrieved refractive in- 
dex averaged over the two modes is consistent with the 
AERONET refractive index. The retrieval results of RSP 
measurements with different spectral ranges and angular res- 
olutions confirm the findings of the synthetic study, i. e., that 
five viewing angles is sufficient for aerosol retrievals and that 
the inclusion of measurements at 1590 nm is important. 
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Appendix: Models for surface reflection of radiance and 
polarized radiance 


The Rahman-Pinty-Verstraete model is used with the fol- 
lowing form: 


(COS !?0 COS & Y ) k 1 

R I (k,#v,#o,4>) = 7 — ^ ^-n-poW F (y) 

(COSJ>0 +COS!7y) 

(l + tf(G)), 

1 — g 2 

F ^ (1 + g 2 -2gcosy)*- 5 ’ 

, , , , 1-Po 


1 + R(G) = 1 + 


1 + G’ 


G = J tan 2 #o + tan 2 i? v — 2 tan #o tan | $ v | cos 0. 


Here po, 8 and k are three independent parameters; X is wave- 
length; y is scattering angle; 0 is the difference of azimuth 
angles, 0 = 0 V — 0o; # v and t?o are viewing and solar zenith 
angles, respectively; F(y) is Henyey-Greenstein phase func- 
tion; 1 + R(G ) is an approximation of the hot spot effect. 


The modified Fresnel model is written as 


an F 0 (m, y) 

R v (& v ,&o,(p) = — 7 — — /(u v ,«o)/sh(p), 

4 U n (COS T7 V + COS 17 o) 

/ 1 +cosA: K (7r — y)\ 3 
/sh(y) “ ( 2 ) 


( 5 ) 

(6) 
( 7 ) 


Here /(n v ,«o) is a Gaussian function which describes 
the distribution of facets over orientation, where a is 
the standard deviation of surface facet slopes. / s h(y) is 
a shadowing function with free parameter k y . F p (m,y) 
is the element F 21 of the Fresnel scattering matrix with 
refractive index m. a is a scaling parameter, no and 
n v are vectors with form (sin#ocos0o, sin#osin0o, cos0o) 
and (sin |# v | cos0 v , sin |z? v | sin0 v , cos0 v ). u n is defined by 

|w V +« 0 | ’ w here n\ and n q are z components of the vector n v 
and iiq. 
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